Optimal control of the silicon-based donor electron spin quantum computing 
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We demonstrate how gradient ascent pulse engineering optimal control methods can be imple- 
mented on donor electron spin qubits in Si semiconductors with an architecture complementary 
to the original Kane's proposal. We focus on the high-fidelity controUed-NOT (CNOT) gate and 
explicitly find its digitized control sequences by optimizing its fidelity over the external controls of 
the hyperfine A and exchange J interactions. This high-fidelity CNOT gate has an error of about 
10~^, below the error threshold required for fault-tolerant quantum computation, and its operation 
time of 100ns is about 3 times faster than 297ns of the proposed global control scheme. It also re- 
laxes significantly the stringent distance constraint of two neighboring donor atoms of 10 ~ 20nm as 
reported in the original Kane's proposal to about 30nm in which surface A and J gates may be built 
with current fabrication technology. The effects of the control voltage fluctuations, the dipole-dipole 
interaction and the electron spin decoherence on the CNOT gate fidelity are also discussed. 

PACS numbers; 03.67.Lx, 82.56.Jn, 85.35.Gv 



One of the important criteria for physical implemen- 
tation of a practical quantum computer is to have a uni- 
versal set of quantum gates with operation times much 
faster than the relevant decoherence time of the quan- 
tum computer. In addition, high-fidelity quantum gates 
to meet the error threshold of about 10"'* (recently shown 
to be about 10""^ in [l|) are also desired for fault-tolerant 
quantum computation (FTQC). There have been several 
different approaches in optimal control of quantum gate 
operation problems 0, H, 01 ■ This work focuses on find- 
ing control parameter sequence in near time-optimal way 
using the gradient ascent pulse engineering (GRAPE) Q 
approach for a high-fidelity CNOT gate in Si:P based 
donor spin quantum computer architectures [EJl, S [1| 
where the electron spin is defined as qubit 0. The 
GRAPE approach partitions a given time into several 
equal time steps, and in each time step of the sequence, 
the amplitudes of control parameters are set to be con- 
stant. For a desired operation, we can define the trace 
fidelity between the desired operation and the unitary 
operation from the sequence. Since we can calculate the 
derivative of fidelity with respect to the control ampli- 
tudes (gradient ascent) in each step, we will be able to 
obtain, given the required fidelity, the near time-optimal 
control sequence numerically. Recently, the GRAPE al- 
gorithm has been applied to the coupled Josephson qubit 
quantum computing and the numerically optimal 
control time for a CNOT gate is found to be 55ps [1| 
instead of 255ps in Ref. [Toj . 

The architecture of Si-based donor spin quantum com- 
puter 0, 0, [1| is composed of '^^P atoms doped in a 
purified ^®Si host where each phosphorus has an electron 
spin and a nuclear spin. In a constant magnetic field Bq 



applied in the z direction, the single-qubit Hamiltonian 
can be written as H = gefJ,BBocr'^/2—gn^nBoa^^/2+Acr'^- 
cr", where the effective electron g- factor in Si = 2, the 
g-factor for a ^^P nuclear spin gn = 2.26, and the hyper- 
fine interaction A « 1.21 x 10~''eV. According to numeri- 



cal calculations llj , it may be possible to vary the hyper- 
fine interaction with A-gate voltage by up to « 50% be- 
fore the donor electron is ionized. Similar to the globally 
controlled electron spin quantum computing scheme [9|, 
we apply a microwave (MW) magnetic field Bac to allow 
for X-axis rotations and also always keep the Bac field on 
as it may not be easy to control and turn on/off the Bac 
field quickly at the precise times in experiments. If we ini- 
tialize the nuclear spins to the spin up state p^ . we can 
use the energy states of |TeTn) and lietn) as a qubit [§]. 
Following Ref. [1], by defining uj{A) = AE{A)/h, where 
A^(^) = 5e/iB5o + 2^-t-[2^V(5eMBSo/2 + 5„^„Bo/2)], 
we obtain the reduced Hamiltonian in the frame rotating 
with the MW field 



H = hAu 0-2/2 + geHBBac(Jx/2, 



(1) 



where Alu = uj{A)- 



and LOac is the angular frequency 



of the MW field Bac- We tune ojac to be the electron 
spin resonance frequency obtained when no voltage is 
applied to the corresponding A gate, i.e., Uac = 
Then the qubits will effectively rotate around the x-axis 
when Auj = (or equivalently A ~ Aq), and around an 
axis which is slightly tilted when Aoj 7^ (or A Aq) 
described by Eq. ([T]). 

The effective reduced two-qubit Hamiltonian, approx- 
imated from assuming that the nuclear spins are frozen 
out to be always up, in the rotating frame is then 



H= hALUial/2- 



hAuj2ol/2- 



gcliBBac{(yl 
— 2e 



^)/2 



(2) 
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where J is the exchange interaction between two adjacent 
donor electron spins. We will use the reduced Hamilto- 
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FIG. 1: (Color online) Fidelity versus time for the CNOT 
gate, (a) gives the trace fidelity against time, while (b) shows 
deviation logj^Q(l — Ftr) from fidelity. 

nian to obtain control sequences by optimizing the fi- 
delity of CNOT gate operations using the GRAPE ap- 
proach. Simulations on the full two-qubit Hamiltonian, 

+geflBBacCOaUJacti(Tl'' + <7l'')/2 + AiO-^" ■ CT^" 

+gefiBBacSmujactial'' + af)/2 + A2a^' ■ (t^" 

-gnf^uBac sinojacticrl'' + uf)l2 , (3) 

with the control sequences found will also be performed 
for error comparison. 

Since the Bac field is always on in this scheme, elec- 
trons will undergo a rotation around the x-axis when 
there are no voltages applied on A gates, i.e. Aoj = 0, 
with an angular frequency of fio — geP'-BBac/f^- While 
the target electrons perform a particular unitary oper- 
ation within time t, every spectator qubit will rotate 
around the x-axis with an angle of 6^ = ^^t. If 6^ 
does not equal to 2n7r where n is integral, another cor- 
rection step will be required for the spectator qubits. 
Thus it will be convenient to choose the operation time, 
t — 2mT/Qo = 2nHTT/{gefJ-BBac), such that there is no 
need for correction for spectator qubits. The Bac field 
is usually very small compared with the Bq field. For a 
given time t, we choose rt = 1 in the reduced and full 
Hamiltonian simulations. In this case, when the con- 
trol duration is 100ns and n = 1, the strength of Bac is 
3.56 X IQ-^T. 

We first try different piecewise constant control steps 
and numerically calculate in the GRAPE approach the 
fidelity (error) against the time needed to implement a 
CNOT gate with stopping criteria of error in the opti- 
mizer set to 10~^ in order to economize the simulation 
time. Here, the error is defined as 1 — Ftr, where Ftr 
is the trace fidelity defined as Ftr = |Tr{J7|)J7F}P with 
Ud being the desired unitary operator in a given time t, 
and Up being the optimal unitary operator constructed 
by our control sequence. For each trying value of time t, 
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FIG. 2: (Color online) Near time-optimal CNOT gate control 
sequence with 30 steps in 100ns obtained using the reduced 
Hamiltonian. In (a) and (b), the maximum energy difi'er- 
ence of term from detuning the hyperfine interaction is 
{l/2)Auj/2-K = -14.7MIIz. In (c), the maximum electron- 
electron exchange energy is J/h — 19.96MHz. 



we divide the sequence into 30 piecewise steps, starting 
with each of the initial control amplitudes (^i, A2 and 
J gates; or equivalently Auji, Alu2 and J) by assigning 
a random value to every five steps in time and using a 
cubic spline to fill in the amplitudes of the intermediate 
time steps. The values of the control amplitudes Ai and 
A2 are varied between Ao/2 and Aq [l,|Tlj, and the value 
of J is varied between and Jq, where Jq is chosen for 
the donor separation to be around 30nm. The fidelity 
against time obtained from the optimization of the re- 
duced Hamiltonian Eq. ^ is shown in Fig. [1] In Fig. [1] 
(b), the error is less than 10~^ for times longer than 
100ns, and it is found that 30 piecewise constant control 
steps for the CNOT gate operation will be sufficient to 
meet the required fidelity (error) and the performance 
would not be improved further with more steps. With 
the operation time t — 100ns and stopping criteria of er- 
ror set to 10~^^, we can find that the near time-optimal, 
high-fidelity CNOT gate control sequence has an error of 
1.11 X 10^^^. The digitized sequence of controls is shown 
in Fig. [21 In a typical Kane quantum computer's scheme, 
the typical value oi J/h k, 10. 2 GHz, which requires the 
separation between two neighboring donors to be about 
10 ~ 20nm [^J. This sets a stringent fabrication condition 
to fabricate surface A and J gates within such a short dis- 
tance. One of the great advantages in our scheme is that 
the maximum exchange energy in our simulation is only 
J/h « 20MHz. This corresponds to a donor separation 
around 30nm [l^. To fabricate gates of this size is 
within reach of current fabrication technology. 

We next apply the control sequence of the CNOT gate, 
obtained from the optimization of the reduced Hamilto- 
nian Eq. ([21), to the full spin Hamiltonian Eq. We 
simulate the CNOT gate numerically with initial four 
different computational basis electron spin states, |00)g. 
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FIG. 3: (Color online) Time evolution of the CNOT gate in 
the rotating frame, simulated using the full Hamiltonian with 
4 different initial electron-spin input states. All the nuclear 
spins are initially spin-up. 



TABLE I: Summary of the CNOT gate errors. 



Input state, 


Expected output 


Error 


Probability that 




5 |00>„ 


state, ® |00)„ 


(1 -p)" 


nuclear spins flip' 


|00>,(5 


5 |oo)„ 


|00),®|00)„ 


1.80 X 10^** 


1.57 X lO"'^ 


|01>e€ 


5 |oo)„ 


|01>e®|00)„ 


1.80 X 10"^ 


2.00 X 10"^ 




5 |oo)„ 


|ll>e®|00)„ 


1.92 X lO"*^ 


1.93 X 10"*^ 


|ll>e€ 


5 |oo)„ 


|10>,®|00)„ 


1.20 X 10"^ 


1.56 X lO"** 



"The output reduced density matrix of the electron spins is ob- 
tained by tracing over all the nuclear states. 

'Here, we trace the total output density matrix over the electron 
spin states to obtain the reduced density matrix for the nuclear 
spin states to compute the flipping probability. 



dies off as where d is the distance between two 

qubits. The dipole-dipole interaction Hamiltonian can 
be written as 



Hd = D [ct^" ■ a^" ~ 3(cri'= . ^)(^2e . 



(4) 



|01)g, |10)g, and but the same nuclear spin-up 

state, where \0)^ means the electron spin is up. The 
final reduced electron density matrix is defined as the 
composite density matrix traced over all the nuclear spin 
states. The errors of the full Hamiltonian CNOT gate 
operations with the four input electron spin basis states 
evolving to their correspondingly expected output elec- 
tron spin states are shown in Table HI Here the error 
is defined as 1 — P, where P is the probability that the 
qubits are in our desired quantum state after the CNOT 
operation. The time evolutions of the states of the CNOT 
gate are shown in Fig. [31 The error are about 10~^ which 
are below the error threshold 10"'* (10"'^ in [ij) required 
for FTQC. Most of the errors result from the accuracy of 
the second-order approximation in A of Eq. ^ since the 
hyperfine interaction A would cause both electron spins 
and the nuclear spins to flip in the full Hamiltonian 
The CNOT gate operation time of 100ns is about 3 times 
faster than the globally controlled electron spin scheme 
[il of 297ns [in g, the indicated CNOT time is 148ns 
that is due to a factor of 2 missing in the denominator of 
their Hamiltonian Eq. (6)]. The error probabilities that 
nuclear spins may flip after the CNOT gate operation for 
the four input electron states are around 10~^ (see Table 
H]). If we repeat the CNOT process N times by simply 
inputting the same pure electron state \ij) but not reini- 
tializing the nuclear state each time, the errors of the 
CNOT gate operations will accumulate. The numerical 
results indicate that in the worst case of the electron spin 
input state |10)g, after around 60 (250) times of opera- 
tions, the error sums up to 1.03 x 10""* (0.79 x 10~^). 
Therefore in order to maintain FTQC, one has to reini- 
tialize the nuclear spin state before about 60 (250) times 
of operations. 

Although the exchange interaction dies off exponen- 
tially with distance, the dipole-dipole interaction that 
couples every pair of electronic spins in the system only 



where D — '^iQ^^a is the dipolar interaction energy, 
7e — is the gyromagnetic ratio of the electrons, 

and n is the unit vector in the direction joining the 
two electrons. In our scheme, the separation of the 
two donor qubits is around 30nm, and thus the corre- 
sponding D K, 1.98 X 10~^^eV, which is still flve or- 
ders of magnitude smaller than the exchange energy J 
used in our scheme. We simulate the optimal control 
sequences obtained previously with the full Hamiltonian 
plus the dipole-dipole interaction Hamiltonian to see its 
effect. Since the first term in Eq. ([4]) has the same 
form as the exchange energy, we may combine this term 
with exchange energy. So what we really need to care 
about is only the second term of Eq. ^ , which becomes 



Gy with the donors aligning along 



the h = y axis. The fidelities of the simulation results 
are slightly worse than the case without dipole-dipole in- 
teraction, but they are almost the same and the errors 
are still below the error threshold lO""' (lO^'^ in [l|) re- 
quired for FTQC. So the dipole-dipole interaction may 
dominate for larger separations, but it is still too small 
to decrease significantly the fidelity of the CNOT gate 
operation. 

Since we apply voltages on the A and J gates to con- 
trol the strengths of hyperfine interaction and exchange 
interaction, there might be noise induced from the (ther- 
mal) fluctuations in the control circuits, which then cause 
the uncertainties of the control parameters and decrease 
the fidelity of a specific operation. To analyze the de- 
crease of fidelity due to these uncertainties, we model the 
noise on the control parameters Ai, A2 and J as inde- 
pendent white noise with Hamiltonian written as — 

TAUt>^^ ■ + ^Ai2{tW' ■ + TjUt)'^^'' ■ 

where the mean of the continuous time random pro- 
cesses {£,i{t)) = 0, the correlation functions {£,i{t)^j{t')) ~ 
6ijd{t—t'), and F^ and Fj are the spectral densities of the 
noise signals, which have the dimension of (energy)2/Hz. 
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FIG. 4: (Color online) Contour plot of logarithmic errors sim- 
ulated under different spectral densities and Fj of the 
white noise signals on the control amplitudes of A and J of 
the full Hamiltonian. The unit of F^/h^ and T^/h^ in the 
plot is Hz and both of the axes are also in logarithmic scales. 

We simulate the optimal control sequence in the pres- 
ence of the white noise through the effective master equa- 
tion approach The contour plot of the logarithmic 
errors of the full-Hamiltonian simulation results due to 
the white noise is shown in Fig. 3] To satisfy the error 
threshold 10"'* (IQ-^ in [H) of FTQC, the spectral den- 
sities, Fj/ft,^ and r\/h^, have to be smaller than 6.2Hz 
and 13Hz (63Hz and 125Hz), respectively. This precision 
of control should be achievable with modern electronic 
voltage controller devices. For example, it was stated in 
that the spectral density of the gate voltage fluctua- 
tions for good room temperature electronics is of order of 
10~i*V'^/Hz, comparable to the room temperature John- 
son noise of a 50 resistor. At a particular bias voltage, 
the gates have a frequency tuning parameter a = df/dV, 
estimated to be 10 ~ 100 MHz/V Therefore, the 
spectral density of energy fluctuations of the control pa- 
rameters for good room temperature devices can be es- 
timated to be lO^"' ^ lO^^Hz that is still much smaller 
than 6Hz required by the error threshold of 10"'*. 

The decoherence time T2 for P donor electron spin in 
purified Si has been indicated experimentally [l5| to be 



potentially considerably longer than 60 ms at 4K. It has 
been shown [8|] that the two-qubit gate fidelity of Kane's 
quantum computer is limited primarily by the electron 
decoherence time, e.g., a typical error of CNOT is 8.3 x 
10~^ with operation time of 16/iS for a simple dephasing 
model of r| = 60ms. In our scheme, the CNOT gate time 
is much faster and we expect the decoherence effect may 
decrease the fidelity less. The error with decoherence can 
be estimated to be 1 — F^e"*/-^^, where Fr and t are the 
trace fidelity and operation time of the gate, respectively. 
For this simple estimate, the error is about 2.7 x 10~^, 
below the error threshold of 10"'' (10"^ in [l|). 

In summary, we have applied the GRAPE approach to 
find the near time-optimal, high-fidelity CNOT gate con- 
trol sequence. A great advantage of the CNOT gate se- 
quence is that the maximum value of the exchange inter- 
action is J/h w 20MIIz which is about 500 times smaller 
than the typical value of 10. 2 GHz in and yet 

the CNOT gate operation time is still about 3 times faster 
than in This small exchange interaction relaxes sig- 
nificantly the stringent distance constraint of two neigh- 
boring donor atoms of about 10 ~ 20nm as reported in 
the original Kane's proposal Q to about 30nm. To fabri- 
cate surface gates within such a distance is within reach 
of current fabrication technology. Unlike traditional de- 
composition method that decomposes general gate op- 
erations into several single-qubit and some interaction 
(two-qubit) operations in series as the CNOT gate in §, 
the GRAPE optimal control approach is in a sense more 
like parallel computing as single-qubit {Ai and A2 both 
on) and two-qubit (J on) operations can be performed si- 
multaneously on the same qubits in parallel (see Fig. [2]). 
As a result, the more complex gate operation it is ap- 
plied, the more time one may save, especially for those 
multiple-qubit gates that may not be simply decomposed 
by using the traditional method. So the GRAPE ap- 
proach may prove useful in implementing quantum gate 
operations in real quantum computing experiments in the 
future. We acknowledge supports from NSC under Grant 
No. 97-2112-M-002-012-MY3, from NTU under Grants 
No. 97R0066-65 and 97R0066-67, and from the NCTS 
focus group program. We are grateful to NCHC for com- 
puter time and facilities. 
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